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Abstract. - There is an ongoing debate about the nature of the bosonic excitations responsible 
for the quasiparticle self energy in high temperature superconductors - are they phonons or spin 
fluctuations? We present a careful analysis of the bosonic excitations as revealed by the 'kink' 
feature at 70 meV in angle resolved photoemission data using Eliashberg theory for a d-wave 
superconductor. Starting from the assumption that nodal quasiparticles are not coupled to the 
(tt, 7t) magnetic resonance, the sharp structure at 70meV can be assigned to phonons. We find 
that not only can we account for the shifts of the kink energy seen on oxygen isotope substitution 
but also get a quantitative estimate of the fraction of the area under the electron-boson spectral 
density which is due to phonons. We conclude that for optimally doped Bi2Sr2CaCu20s-|-<5 phonons 
contribute ~ 10% and non-phononic excitations ~ 90%. 
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Experimental evidence based on the analysis of the op- 
tical properties of the high T c cuprate superconductors 
suggested early on that the charge carriers were strongly 
coupled to a spectrum of bosonic excitations in the 40 
to 60mcV range [1,2]. This was confirmed by the dis- 
covery of a slope change or kink in electronic disper- 
sion of the quasiparticles (QP) by angle resolved pho- 
toemission spectroscopy (ARPES) in the same energy 
range [3-5]. The boson coupling that manifests itself as 
peaks in the real part of the QP self energy extracted 
from the ARPES dispersion curves can be seen in opti- 
cal spectroscopy as features in the optical scattering rate, 
a term in the generalized Drude formula for the optical 
conductivity [1,2,9,12,17,19-21] and in tunneling [27,28]. 
Surprisingly, the nature of the bosons involved remains 
highly controversial [6-21]. The subject is of great inter- 
est since such a coupling, would give us a working model of 
superconductivity in these materials, at least within the 
Migdal-Eliashberg formalism with retarded interactions. 
In support of this approach, recent numerical solutions 
of the Mott-Hubbard model suggest that in the high T c 
cuprates retarded interactions do provide most of the pair- 
ing [22,23], with the energy scale set by the size of the anti- 
ferromagnetic exchange constant J. These calculations do 



not support the suggestion by Anderson [24] that the en- 
ergy scale for the retardation might be much higher and 
set by the Hubbard U which would imply non-retarded 
pairing. Another mechanism considered in the literature 
is kinetic energy pairing for which the kinetic rather than 
the potential energy is lowered as the temperature is low- 
ered below T c . In an Eliashberg formalism this can be 
simulated by a phenomenological reduction in quasipar- 
ticle scattering which could be due to a reduction in the 
electron-boson spectral density [25,26]. But this, presum- 
ably, applies only to the superconducting state and would 
give no dispersion 'kink' above T c . 

The controversy over the bosonic spectrum centers on 
the nature of the bosons: are they phonons or magnetic 
excitations. Neutron spectroscopy does not resolve this 
problem as it shows both phonon and magnetic excita- 
tions in the energy region of interest. What is needed is a 
fingerprint experiment, one that unambiguously identifies 
the suspect bosons. In the past a number of candidate ex- 
periments have been proposed: doping and temperature 
evolution of the bosonic spectrum [5,9,13,14], substitu- 
tion with magnetic and non magnetic impurities [15], but 
the gold standard for the resolution of this problem is the 
oxygen isotope effect. A spectral feature caused by a lat- 
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tice vibration involving oxygen will shift its frequency in a 
typical oxide by 6.5% on 16 — > 18 substitution whereas 
a magnetic excitation will yield little or no change [29] . 

There have been many studies of the effect of oxygen iso- 
tope substitutions on the boson structure in the cuprates. 
While the early ARPES experiments could not resolve the 
expected 4.4 me V frequency shift in the kink, an optical 
study by Wang et al. [30] on the frequency of the shoulder 
in the reflectance of undcrdoped YBa2Cu307_<5 found no 
evidence of an isotope effect. Recent scanning tunneling 
microscope (STM) work has shown a shift in the 52 meV 
structure of about the expected amount for an oxygen 
mode [27], but this may be explained in terms of inelastic 
tunneling through a barrier [31-33]. Very recently new 
high precision low-energy ARPES data of Iwasawa et al. 
[34] on a high quality single crystal of optimally doped 
Bi 2 Sr 2 CaCu 2 8+( 5 (Bi2212) with a T c of 92 K have re- 
vealed a shift of the kink at 69 meV seen in the dispersion 
curves in the nodal direction for 16 O — > 18 O substitution 
with little change in the spectrum outside the immedi- 
ate kink region. By providing a more detailed analysis of 
the spectra in Ref. [34] we will show here that not only 
do these new measurements provide a fingerprint for the 
presence of a phonon contribution to the self energy of the 
charge carriers, but that they permit us to calculate the 
fraction of the bosonic coupling parameter A associated 
with phonons and the fraction with magnetic excitations 
and conclude that the phonon contribution is ~ 10% of the 
area under the spectral density. The major contribution 
to the glue is non phononic. 

A key to an understanding of the bosonic structure in 
ARPES is the recognition that there are several features 
that can contribute to the kink in the dispersion curves 
and that the challenge is to separate their contributions. 
We will show that the new experiments of Iwasawa et 
al. [34] allow us to do that unambiguously. Calculations 
of scattering by spin fluctuations suggest that an impor- 
tant role is played by bosons that have a (tt, tt) wave vector 
and they are particularly efficient in scattering charge car- 
riers near the antinodes but have little effect on the nodal 
quasiparticles. Neutron scattering experiments show that 
the local (q averaged) susceptibility has a strong peak in 
energy around 40 meV, but also that the q dependent sus- 
ceptibility is peaked at (tt, tt) in this energy range. As a 
result, in frequency dependent measurements, the interac- 
tion with this magnetic peak is small in the nodal direction 
because Fermi surface to Fermi surface electronic transi- 
tions involving the spin one resonance with momentum 
transfer of (tt, tt) cannot take place due to the particular 
Fermi surface geometry but this becomes more possible 
as we move towards the antinodal direction. This agrees 
with ARPES measurements that show the growth of a 
sharp feature as the momentum of the quasiparticles ap- 
proaches the zone boundary at (0,7r) [10]. 

What then scatters the quasiparticles in the nodal direc- 
tion? We suggest that spin excitations that are not cen- 
tered at (tt, tt) but spread out as a continuum and phonons 



that are also expected to be scattered more or less inde- 
pendent of QP vector q. ARPES measurements [18] show 
that while there is a sharp spectral feature in the 65 meV 
region for nodal quasiparticles it is weak. Nevertheless, as 
we have argued above, it cannot be assigned to the (tt, tt) 
spin resonance and we suggest this feature to be the resid- 
ual phonon scattering superimposed on a broad magnetic 
background. In what follows we will show that this idea 
is fully consistent with all the data of Iwasawa et al. [34] 
and thus provides us with a quantitative tool to estimate 
the relative contribution of phonons and spin fluctuations 
to high temperature superconductivity. 

The changes in the real and imaginary parts of the 
self energy with isotopic substitution are best illustrated 
with a simple model for the effect of a single Einstein 
mode at frequency lue on the self energy. To simplify 
the algebra we first show the results for a mode of zero 
width. The electron-boson interaction spectral density 
I 2 x(u) = A5(lu — cue) where A is the area under the Dirac 
5-function 8(x). The resulting QP self energy in the nor- 
mal state at zero temperature (T = 0) is 



T,(uj) = Si(w) + iS 2 (w) = Aln 
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where 9(x) is a step function. The dotted and the dash- 
dotted curves in Fig. 1(b) illustrate these functions where 
we have modified the spectral density by giving it a 
Lorcntzian width. We note that under isotopic substi- 



O, the frequency of the mode is reduced 



tution 16 O -> 18 
to 7^b, but the coupling constant too is reduced 

to A -> 7^4 = ^ iso where 7 = ^16/18 = 0.94. The mass 
enhancement parameter A = 2A/uje remains unchanged. 
In the limit lu ^ lue the real part of the self energy is 
linear in to with slope A. As Fig. 1 shows, isotopic substi- 
tution leaves two 'fingerprints' in the self energy spectra, 
the well known shift of the mode frequency by 6% but a 
second, less known one, a reduction in amplitude of 6%. 
In the real part of the self energy, Fig. 1(a), the peak is 
shifted to lower energy and reduced in amplitude by the 
same factor while for the imaginary part, Fig. 1(b), the 
onset of scattering is also moved to lower energies and its 
saturation value is reduced by the same 7. The inflection 
points are found at lol and LdL.iso = ju>l respectively. 

Many of the features of the model of a single bosonic 
mode shown in Fig. 1 do not agree with recent data of 
Iwasawa et al. [34] even when we account for the changes 
in QP properties that are brought about by superconduc- 
tivity. First, coupling, predominantly to a single mode 
peaked around lol, whatever its origin, leads to a satura- 
tion in the magnitude of the QP scattering rate, T, 2 (uj), at 
frequencies above lu^. Such a saturation is not seen in the 
published ARPES data obtained by fitting a Lorcntzian 
form to momentum distribution curves, nor is it seen in 
optical data. Both methods find scattering rates [16, 17] 
which continue to increase up to energies of order 400 meV. 
This is well beyond any phonon energy (~ 86 meV in 
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Fig. 1: (Colour online) (a) Minus the real part of the QP 
self energy, — £i(w) ( ie O solid and ls O dashed lines) vs en- 
ergy uj. The two thin dashed and dotted vertical lines in- 
dicate the position of the peak in the I 2 x(&) spectrum, 
and u)-^ - 1SQ , respectively, (b) The same as (a) but now for 

—£2(0;). Lorentzian model I 2 x{u) spectrum for ls O (dashed- 
dotted line) and 7 2 Xi so (<^) for 18 (dotted line) right hand 
scale. 



Bi2212 [35]) and shows that in addition to the sharp 
bosonic peak the charge carriers are also coupled to high 
energy boson modes of some other origin. Secondly, the 
model predicts a reduction of the amplitude of the real 
part of the self energy of 6% but the data of Iwasawa et al. 
shows little change in the range around 70 meV. While the 
frequency of the peak is shifted by the expected amount, 
the magnitude of the maximum in — Ei(w) appears not 
to be significantly changed. This contradicts the single 
phonon model of the self energy. 

We will next show that these contradictions can be re- 
solved by assuming that in addition to a (weak) phonon 
contribution there is a strong contribution from magnetic 
scattering and that there is additional influence due to fi- 
nite band widths. Taking this into account we are able to 
estimate the relative strength of the phonon and magnetic 
contributions to the QP self energy. 

To accomplish this we use the recent high resolution 
ARPES work of Zhang et al. [18]. These data were an- 
alyzed using maximum entropy inversion [19] in the su- 
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Fig. 2: (Colour online) (a) Presents hts to the real part of the 
superconducting state QP self energy ARPES data of Ref. [18] 
as reported by Schachinger and Carbotte [19] for T — 17 K 
(solid line and open diamonds) and 45 K (dashed line and open 
squares), (b) The I 2 x(^) spectrum from maximum entropy 
inversion of the T — 17 K Zhang et al. [18] data (solid line) 
with A = 1.19, and the area under the spectrum A = 51.9 meV. 
The (red) dashed line shows the spectrum I 2 xi so (^) which was 
used to simulate the isotope effect (depicted also in the inset 
for clarity). Only the area of the peak equal to 6.15 meV or a 
partial A of 0.21 is shifted in energy by 6% with A = 1.19 and 
A iso = 51.4 meV. 



perconducting state based on finite band width c?-wave 
Eliashberg equations, in order to determine the spectral 
density I 2 x(u) of the bosons contributing to the nodal 
direction QP self energy. In Fig. 2(a) we show our fits, 
solid curve for T = 17 K (dashed curve for T = 45 K) to 
the ARPES data, open circles (open squares). Results for 
the electron-boson spectral density I 2 x(uj) at T = 17 K 
and a bandwidth of 1.2 eV are shown in Fig. 2(b) as the 
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Fig. 3: (Colour online) The real part of the QP self energy 
— Ei(a>) vs energy ui of Bi2212 in the superconducting state at 
T — 17 K. The solid line corresponds to the solid line in the 
inset of Fig. 2. The dashed line is the result due to the isotope 
effect simulated by the spectral function I 2 X[ S0 (u)- The peak 
position shifts down by 4.5% and its amplitude by 3%. (b) The 
same as (a) but for a cutoff of 200 meV in the spectra of Fig. 2. 
The dashed-dotted and dotted lines correspond to the 16 O and 
ls O data for cut zero presented by Iwasawa et al. [34] in their 
Fig. 2(b) rescaled to meet our theoretical results in amplitude. 

solid (black) line. This J 2 x(w) spectrum is very differ- 
ent from the single sharp Lorentzian form used so far. It 
has a broad peak around uj p = 65 meV superimposed on 
a large background extending to 400 meV. The cutoff in 
this spectrum depends critically on the choice of the bare 
dispersion curve made in the ARPES study. The renor- 
malizations end at the crossing between bare and dressed 
dispersions and in the work of Zhang et al. [18] this was 
taken to be 400 meV. This choice is consistent with optical 
data [12.17,20,21]. It is also consistent with the obser- 
vation that the QP as well as optical scattering rates still 
increase with increasing u> in the range of a few 100 meV. 

Our model of the isotope substitution on the real part of 
the self energy is based on shifting only the peak in I 2 x(uj) 
shown by the (red) dashed curve in Fig. 2(b) leaving the 
background unchanged as depicted in the inset, and re- 
sults for the self energy are shown in Fig. 3(a). The solid 
curve is for 16 O and the dashed curve for 18 O. The results 



are in good overall agreement with the experimental data 
of Iwasawa et al. [34]. For instance, we predict a shift in 
the position of the peak of the order of 4.5% against the 
observation of 5 ± 0.8 %. We also predict a reduction in 
the amplitude of the self energy maximum of ~ 3% some- 
what larger than is observed in the data. Note that the 
position of the peak in — Si(oj) is not trivially correlated 
to the position of the peak in the spectrum, lo v (ujp.iso), 
nor to their values displaced by the gap amplitude, uj p + A 

Figure 3(b) presents the results of a numerical simula- 
tion using a smaller energy cutoff of 200 meV as done by 
Iwasawa et al. [34]. The solid and dashed lines are for 
16 and 18 0, respectively. The temperature is T = 17 K 
and the bandwidth had to be reduced to 0.6 eV. We also 
added experimental data for 18 O (dashed-dotted line) and 
for 16 O (dotted line) from Fig. 2(b) (cut zero) of Iwasawa 
et al. rescaled to meet our Ei(w) data in amplitude. The 
qualitative agreement between data and numerical simu- 
lation is now even better. As a result of the new cutoff 
the peak position in our theoretical results moved down 
by ~ 6meV. Also the difference in amplitude is now less 
than one meV and this explains why almost no difference 
in amplitude can be seen in the data. As a result of this 
simulation we can safely conclude that the data can be re- 
produced by a 7 2 x(w) spectrum which has the same shape 
as the spectrum presented by the solid line in Fig. 2 with 
the broad peak which will be shifted by isotope substitu- 
tion now centered between 55 and 58 meV. 

The area under the electron-boson spectral density 
which we have shifted in energy by 6% corresponds to 
6.13 meV or about 10% of the total area under the I 2 x{^)- 
If we attribute this shift in area entirely to an electron- 
phonon coupling, the mass enhancement factor involved 
is A ~ 0.2. This value is much smaller than estimated in 
Ref. [34] but is of the order found in band structure calcu- 
lations on related cuprates [36-38] . It also compares favor- 
ably with the total value of mass enhancement X to t = 0.23 
found in studies [39] of anisotropic electron-phonon cou- 
pling due to the oxygen buckling mode (~ 36meV). The 
coupling to the breathing mode (~ 70 meV) was found to 
be much smaller: X tot ~ 0.02. Our value of A = 0.2 corre- 
sponds to less than 20% of the total A found through max- 
imum entropy inversion [19] of the nodal direction ARPES 
data of Ref. [18]. 

In summary we analyzed recent nodal direction ARPES 
data of the effect of 16 O — ► 18 O substitutions on the QP 
self energy. We find that the data are consistent with a 
model where only ~ 10% of the self energy is through 
coupling to lattice vibrations, so that the major part of 
the glue function is of another origin, we believe spin fluc- 
tuations. We conclude that the new measurements need 
not imply that phonons play a large role in the super- 
conductivity of these materials. If the mass enhancement 
parameter A is used instead of the area under I 2 \(ui) as a 
measure of the strength of the electronic renormalizations, 
the phonon contribution to A is 0.2 or less than 20% of the 
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total A = 1.19 for our realistic model for the spectral den- 
sity in Bi2Sr2CaCu208+,5 obtained from inversion of nodal 
direction ARPES data. 
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